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ABSTRACT 

We discuss dissipative processes in st rongly Kyrotropic , nearly collisionless plasma in clusters of 
galaxies (ICM). First, we point out that Bra ginskvl (|1965| ) theory, which assumes that collisions are 
more frequent that the system's dynamical time scale, is inapplicable to fast, sub- viscous ICM motion. 
Most importantly, the electron contribution to coUisional magneto-viscosity dominates over that of 

/ \ 1/4 

ions for short-scale Alfvcnic motions with wave length satisfying I < (^j ~ Ikpc (where A 

is particle's mean free path, /? is the plasma pressure parameter and TOe,p sue electron and proton 
masses). Thus, if a turbulent cascade develops in the ICM and propagates down to scales < 1 kpc, it 
is damped coUisionally not on ions, but on electrons. 

Second, in high beta plasma of ICM, small variations of the magnetic field strength, of relative 
value ~ 1//3, lead to development of anisotropic pressure instabilities (firehose, mirror and cyclotron). 
Unstable wave modes may provide additional resonant scattering of particles, effectively keeping the 
plasma in a state of marginal stability. We show that in this case the dissipation rate of a laminar, 
subsonic, incompressible flows scales as inverse of plasma beta parameter. We discuss application to 
the problem of ICM heating. 
Subject headings: galaxies: clusters: general 



1. INTRODUCTION 

One of the key problems in physics of interclus- 
ter medium (ICM) is the absence of strong cooling 
flows at the centers of galaxy clusters (see, e.g. , 
IPeterson fc FabianI l2006l for a review). It has been 
proposed that heating of ICM by Active Galactic Nu- 
clei (AGN s) may be sufficient to offset the cooling {e.g. 
lBegelmanl f2004'). While the total energy budget of AGNs 
is, in principal, sufhcient to offset the radiative cooling, 
details of how this is achieved are far from clear. 

The observational confirmation of the AGN heating 
model comes from ubiquitous presence of AGN blown 
bubbles, identified by decreased X-ray e mission in Chan- 
dra and XMM maps (McNamara '2000'). These bubbles 
expand and rise in the cluster potentials transferring part 
of their energy to the internal energy of ICM. It has been 
suggested that this process can be very efficient, so that 
a large fraction of the power released by AGN ends up 
as internal energy of ICM. 

The high efficiency of energy dissipation is far from 
obvious. What is required is a distributed increase 
of the entropy of the gas, not jus t of the internal or 
bulk energy (entropy floor problem iLlovd-Davies et aI7\ 
[20001) . The main problem is that these AGNs blown 
bubbles expand, typically, subsonically, as is indicated 
by the general absence of shock signatures ahead of the 
bubbles. In laminar flows at small Reynolds numbers. 
Re < Rccrit ~ 10— 100, dissipation efficiency is oc 1/Re. 
For Re ~ 100, such a low efficiency puts unreasonable 
demands on AGN luminosity. 

2. COLLISIONAL DISSIPATION IN GYROTROPIC 
PLASMA: GYRORELAXATIONAL HEATING 

Ion Larmor radii in ICM, ^ 10^ cm, is some fif- 
teen orders of magnitude smaller than the system size, 
L ^ hundreds of kpc, and Coulomb mean free path. 



A ^ 10 — 30 kpc, for a typical density n 10 ^ cc, 
magnetic fi elds ~ 1 — 10/iG and te mperatures in the keV 
range {e.g. ICariUi fc Tavloij|2002D . Thus ICM is weakly 
collisional, «C A. In addition, it is weakly magnetized, 
in a sense that magnetic fields energy is smaller than 
plasma pressure, (3 — ^ttP/B"^ > 1. Below we will refer 
to this regime as a strongly gyrotropic plasma. 

Dissipation in a strongly gyrotropic plasma proceeds 
in a qualitatively different way from the isotropic case, 
as is exemplified by so called gyrorelaxational heating. 
If in an initially pressure-isotropic plasma the absolute 
value of magnetic field oscillates with a frequency uo and 
relative amplitude 5 — 5B / Bq, then the dissipation rate 
a (so that energy of a parti cle £ changes a ccording to 
d£/dt = a£) in a cycle is feg lBorovs 

a « -— (1) 
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where Vc is collision frequency. Dissipation of energy 
occurs both due to electron and ion collisions, so that 
a = a,, + cti, calculated with corresponding collision 
frequencies v^. and Vi. In the high collision frequency 
regime, i^c 2> tiJ, Eq. (HI) approximates Braginsky's 
result, a oc l/i>c Braginsky (1965). Since ions have 
smaller collision frequency, dissipation in this limit is 
dominated by ions. On the other hand, for rare colli- 
sions, ^, dissipation rate is proportional to collision 
frequency, a (X v^, and is thus dominated by electrons for 
iew > |(me/mi)^/4 K. 0.2 

Consider sub-viscous turbulent motion of ICM occur- 
ring on scale I smaller than mean free path A and medi- 
ated by Alfven waves, so that a typical wave frequency 
is w = Va/1 ^ Cs/{-\fpl). Then for waves satisfying 
teUJ > |(me/mj)^/'*, or for 

l< -=[ — 



1/4 



Ikpc. 



(2) 
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electron viscosity dominates over ion viscosity. For nu- 
merical estimates we assumed T = 10^ K, n — 10~^ cc, 
-B = 5 /i G, so that /3 ~ 10 and mean free path A = 23 
kpc. 

Thus, if a turbulent cascade develops in the ICM 
and propagates down to scales < 1 kpc, it is damped 
coUisional l y not on ions, but on electrons. Thus, 
\Braains'k'i\ Hi 96 A ) theory, which assumes frequent colli- 
sions, tcoUijJ ^ 1, is inapplicable to fast, sub-viscous ICM 
motion. 

3. HEATING IN A BOUND ANISOTROPY MODEL 

Besides binary collisions, plasma can be heated 
through development of electromagnetic turbulence 
which resonantly scatters particles and, thus, provides 
an additional dissipation mechanism. In this Section we 
describe such mechanism of dissipation through develop- 
ment anisotropic plasma instabilities. 

3.1. Viscosity due to binary collisions in a gyrotropic 

plasma 

When Coulomb collision frequency Vc is much smaller 
than cyclotron frequency, lob/vc S> 1, plasma viscos- 
ity is strongly an isotropic, determined by seven coefh- 
cients ()Braginskvl [r965) . In the limit lob ^ oo and slow 
changes of magnetic field, w <C t'c, the only remaining 
coefficient is r?o, responsible for the viscosity along the 
field lines. In this case the viscose stress tensor becomes 
(jlandau fc Lif shits 1982) 



?7o {'ibibj - 6ij) bibkdiVk 



-divv 



(3) 



where bi is a unit vector along the local magnetic field, 
?7o = p/i^c, P = {P\\ + 2Pj_)/3 is total pressure. Below 
we concentrate on the incompressible limit, divv — 0, 
which eliminates reversible compressional heating. For 
incompressible plasma without conductivity, using Eq. 
([3]), the volumetric dissip ation and entropy gener ation 
rates due to viscosity are (jLandau fc Lifshitall975l ) 

Using induction equation, dB/dt = (BV)v, the entropy 
generation rate can be related to the rate o f change of 
magnetic field (jSchekochihin fc Cowlevl[2"006f l: 



dE ^ / I dB 



(5) 



Dissipated power of a gyrotropic fluid is solely due to 
changing magnetic field, which is very different from 
the isotropic case. This result can also be verified if 
we note that in a gyrotropic plasma the entropy is 
S oc (l/2)lnPi_P| (assuming constant density). The 
entropy production is then 



dS 1(P±-P||)2 



dt 



Pi Pii 



(6) 



For binary collisions using P± — P|| = Srjodt In B (Eq 
(fT2l) ). this gives 

^JdtlnBy 



f ^3KlnP)^.o^^ 



consistent with ([5]). 

The differences between the dissipation rates calcu- 
lated using isotropic and anisotropic viscosities can be 
dramatic. For example, for spherical expansion of a bub- 
ble into incompressible fluid, in absence of magnetic fleld, 
the dissipated power is zero (flow field is irrotational) . 
Introduction of a weak (in a sense that /3 3> 1) mag- 
netic field changes this picture completely. In a kine- 
matic approximation (neglecting its dynamical effects, 
so that field lines are just advected with the flow satis- 
fying frozen-in condition) expansion of a bubble into a 
constant magnetic field creates magnetic fields 



Bo = 



sm( 



(1-^-3)1/ 



jj^Bo, Br - - 



COS t 



(l-r')'/'Bo (8) 



where ^ — r/R{t) > 1 and Bg and Br are compo- 
nent of magnetic field in a spherical system of coordi- 
nates aligned with the initial direction of magnetic field. 
Though tangential component of magnetic field diverges 
on the contact ^ = 1 (magnetic draping effect), the 
increase in B-field energy over the initial homogeneous 
field is finite, {1/9)BqR^ and the total heating rate is 
dtE = SrioR^J d^^{dt InBf = 9.5AT]oR^{d\iiR/dt)^. 

This example illustrates an important point: even a 
weak magnetic field may considerably affect plasma dis- 
sipative properties. Inverse situation, when a dissipative 
flow with isotropic viscosity becomes non-dissipative in 
the strongly gyrotropic limit, is also possible. The exam- 
ple is a longitudinal shear, when magnetic field is directed 
along velocity. In the absence of cross-field viscosity there 
is no dissipation. 

3.2. Anisotropic pressure instabilities 

In coUisionless plasma, particles in m agnetic fields 
tend to conserve their adiabatic invariants (jChew et al. I 
fT956h . In case of rare collisions th e equations des cribing 
evolution of pressures becomes (eg lHollwedll985l ) 



dlnP^/B 
Jt 

dlnP||S2 
Jt ^ 



vP\ 
3~ 



2v P\\ - P± 
"3 Pii 



(9) 



(7) 



where P± and P|| are pressure across and along magnetic 
field. 

In a /3 3> 1 plasma, the development of pressure 
anisotropy may lead to firehose, mirror and ion cyclotron 
instabilities when the following conditions are satisfied 

1 

(3\\ — f5j_ > 2, firehose, — > 1 4^ — — mirror (10) 
P|| P_L 

and /3j_//3j| > 1 + fc//3[p^ cyclotron, where (3\\ = SttPw/B"^, 
Pj_ = SttPJB^, 35 < A: < 0.65 and 0.4 < m < 0.42 
iGarv et al. I ()1994D . Cyclotron instability has growth 
rate larger than the mirror instability for /3 < 6 and 
P± > pjj . If initially plasma pressure is isotropic, 
firehose and mirror instability occur when SB/B — 
-2/(3/3o)(firchosc), SB/B = -H/(3/3o)(mirr or) and sim- 
ilar e xpression for ion cyclotron instability (|Garv et al. I 
I1994D : for clarity we do not consider it here. 

The instabilities' increment is maximal at the cyclotron 
frequency, which is very fast compare to any dynamical 
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time. Further change of magnetic field, beyond the limits 
(|3.2|) . will be accompanied by development of instabili- 
ties which will lead to increased scattering rate, either 
due to quasi-linear diffusion or a fully developed turbu- 
lence. As a result, the system dissipates quickly any free 
energy in excess of instability threshold and relaxes to 
the marginally stable state. We expect that the system 
remains at threshold of instability. 

3.2.1. Binary collisions in sub-critical regime 

Binary collisions decrease a level of anisotropy and 
may stabilize plasma. Redefining pressures P± and 
in terms of total pressure p (a trace of the pres- 
sure tensor) and pressure disbalance (Pi — P\\)/p — A, 
= p — |Ap, Pj^ = p + ^Ap, we find 



2pA- 



dB 
dt 



dp 



dA 



dt 
9-3A 



2A2 d\nB 



= 



(11) 



3 dt 

In a /3 S> 1 plasma at the moment of instability A is 
small, |A| ^ 1. For slow changes d/dt <^v this gives 

d\nB , , 

. ... 

This implies that for development of instabilities the dy- 
namical time tdyn ^ l/dtlnB should be relatively short, 
tdynVc < This condition is satisfied by most scales of 
interest in ICM plasma. 

3.3. Dissipation at marginal stability 

As we argued in the previous section, changing mag- 
netic field will lead to development of instabilities that 
will keep the plasma anisotropy at the critical values 
given by Eqs. HU]). Eqs. ^ and the condition pI7)) 
may be regarded as defining an effective scattering rates 

^ " ' ' " ,r = SpdtlnB (13) 



^ef f, firehose — -^PdtlnB, Veff,-, 



At a critically balanced case, the entropy generation 
rate Eq. ^ 



dtS^-dtlnB X 

P 



(14) 



for the firehose and mirror regimes. 

We have arrive at an important result related to effi- 
ciency of dissipation: in a gyrotropic plasma efficiency 
of dissipation is determined not by Reynolds number, but 
by the plasma beta parameter. Typical dissipation time 
scale is (3 times dynamical time, not Re times dynamical 
time. 

The role of effective collisions in energy dissipation in 
a marginally stable regime is, in some sense, opposite 
to the role of binary collisions in a sub-critical regime. 
The entropy production rate and corresponding volu- 
metric dissipated power, Eq. are proportional to 
pressure anisotropy and collision frequency, cx i^(AP)^, 
where AP is the difference in parallel and transverse 
pressures. If the pressure disbalance is due to binary 
collisio ns, then AP oc 1/ v so that the dissipation rate is 
oc l/v ("Bragi nskvlllOBSl ). Thus, before the instabilities 
are reached, increasing collision rate leads to decreasing 
dissipation. On the other hand, for marginally stable 
case AP ~ constant, so that dissipated power is propor- 
tional to the effective collision rate, Eq. 



3.4. Damping of waves at marginal stability 

For Alfvcn waves, perturbations of magnetic field are 
orthogonal to the initial magnetic field, so that variations 
of the absolute value of the field are second order in am- 
plitude. For large enough amplitude, satisfying condition 
[SB/BqY = 6^ > 1//3, this creates conditions favorable 
for mirror instability. The entropy production rate over 
the period is 



dS^ 
dt 



4 



' 2 arccos 



/31 + (52 



(15) 



where the term in parenthesis takes into account phases 
when the amplitude of fluctuations satisfies the mirror 
instability criterion. With Braginsky viscosity, the col- 
lisional damping of Alfven waves is a non-linear effect 
as well, but it has a much steeper dependence of wave 
amplitude and frequency. From ([5]) we find 



dS_ 
dt 



3(5'*w 



4, ,2 



(l + (52)2j,^ 



(16) 



For comparison, in isotropic MHD Alfven waves are 
damped at a rate (|Landau fc Lifshi"t3 ()1982f )) dS/dt = 



4. DISCUSSION 

Our approach follows a long established p rocedure 
of marginal stab ilit y iKennel fc PetschecM (1966); 
Manheimer fc Boris! (|1977D ; iGarv et al. 1 ((19941) : 
Denton et al. I ( 1994J) . when the instability thresh- 
old becomes the limiting value of anisotropy . In 
particular. Quest fc Shapiro ()1996D : iGarv et al. I ()1998D 
applied a bounded anisotropy model to the measure- 
ments of parallel and perpendicular temperatures in the 
solar bow shock region near the Earth magnetosphere. 
It was found that an initial rapid growth of unstable 
waves indeed brings the system back to approximate 
marginal stability. 

What is the relation of the marginal stability condition 
and the conventio nal quasilinear and turbu lence theo- 
ries? According to 'Man heimer fc Boris! (|1977D . both pre- 
dict some level of turbulent fluctuations. Marginal stabil- 
ity approach is applicable if the level of those fluctuations 
is smaller than the one calculated from non- linear theory. 
This, typically, happens when the driver of the instability 
(in our case a large scale motion of ICM plasma) is not 
strong. Assessing whether this is satisfied in case of ICM 
plasma requires full scale calculations of non-linear tur- 
bulence levels, a prohibitively complicated task given the 
uncertainties in both plasma microphysics and details of 
ICM plasma motions. 

The most important effect that was not taken into ac- 
count in the present work is thermal conduction. The 
double-adiabatic equations are valid only when heat fiux 
along magnetic field lines can be neglected. This is the 
main reason w hy the theory rn ay fail {e.g. , the noto- 
rious results of (Kulsrud et al. I [1965'). Neglect of heat 
flux requires that phase velocity of the perturbations 
be much larger than speed of heat carriers, electrons: 
(t) ~ ^Te- This condition may be broken in 

ICM, especially outside of cluster cores. On the other 
hand, enhanced scatterin g rate suppresses conductivity 
(|Levinson fc Eichlei!ll992! ) . The conduction coefficient is 
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K ~ n^v^ ^/vef f (assuming that saturate conductivity 
regime ()Cowie fc McKedll977[ ) is not reached). The ef- 
fective scattering frequency due to development of elec- 
tromagnetic instabilities Eq. (|13p may be higher than 
binary collision rate, so that the conduction coefficient 
will be smaller, k ^ rieV^ / {(3V). Increased scattering 
will also inhibit the onset of saturated regime. 

There is a number of challenges that heating models 
should overcome. Primarily, the heating must be both 
widely distributed and gentle. It is hardly achievable 
with shocks, which provide very concentrated heating at 
the shock location, deposit most of the energy in the 
core and generally contradict the observational absence 
of shock signatures. This, combined with low heat con- 
ductivity in the cores, leads to plasma overheating and 
creation of in yerted entropy gradi ents, contrary to obser- 
vations {e.g. IVoit fc Brvaiill2001[ ). 

The heating in the bounded anisotropy model may be 
distributed. Consider a cluster with a typical density 
profile p oc l/r. Then if bremsstrahlung dominates over 
line emission, the cooling rate is oc (for nearly con- 
stant temperature in the cores). Since an energy flux 



from central source scales as oc as well, this im- 
plies that a heating rate should be independent of a ra- 
dius, and thus independent of the local plasma proper- 
ties. Collisional dissipation clearly cannot produce this. 
On the other hand, if (3 is nearly constant, the heating 
rate will be nearly independent of radius. Thus, at least 
in principle, heating and cooling can be balanced in the 
bound anisotropy model. 

One of the main drawbacks of many simulations of 
ICM is that they use isotropic Spitzer viscosity. Exam- 
ples in ^J3] show that this can produce (at least locally) 
drastically incorrect results, which may either overesti- 
mate or underestimate the real collisional magnetovis- 
cosity (we are not aware of any ICM-relat ed simulations 
with anisotropic viscosity (see, though, Sha rma et al. I 
I2OO6I )). As for the value of the coefficient of viscosity, 
we argued that for binary collision it generally depends 
on electron and ion temperatures and dynamical times 
scales, while in case of marginal stability it is actually 
unrelated to the Spitzer value. Parametrization with re- 
spect to Spitzer may be useful, but we should not put 
too much physical emphasis on it. 
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